Submitted to the Astrophysical Journal 

Preprint typeset using L^T^X style cmulatcapj v. 03/07/07 



O 

(N 

< 

in 



O 

U 

6 
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ABSTRACT 

It is shown that broad-line region (BLR) line profiles ranging from the classic "logarithmic" profile 
to double-peaked, disk-like profiles are readily explained by the distribution of BLR gas proposed by 
Gaskell, Klimek, & Nazarova (GKN) without any need to invoke fundamental differences in the AGNs 
other than differing viewing angles. It is argued that the highly-variable thermal energy generation 
in AGNs originates off axis in regions that cannot be axially symmetric. This off-axis model readily 
explains the varying degrees of temporal correlation found in multi-wavelength variability studies, the 
strong, variable asymmetry of BLR line profiles, the varying time delays in the response of the BLR 
to different continuum events, how narrow velocity ranges of line profiles will often appear to respond 
differently or not at all to continuum variability, complex changes in the Balmer decrement with 
velocity, inconsistent and variable inflow/outflow signatures found in velocity-resolved reverberation 
mapping, the diversity of velocity-dependent polarizations observed, and polarization variability. The 
fundamentally non-axisymmetric nature of AGN continuum variability severely limits what can be 
learned from reverberation mapping. In particular, high-fidelity reverberation mapping is not possible. 
There will be systematic orientation-dependent errors in black hole mass determinations. The effects 
of off-axis emission will mask subtle signatures of possible close supermassive black hole binaries. 
Some tests of the off-axis- variability model are proposed. 

Subject headings: accretion, accretion disks — black hole physics — galaxies: active — galaxies: 
Seyfert — polarization — quasars: emission lines 
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1. INTRODUCTION 



broad-line region (BLR) and 
AGNs have long been enig- 



The nature of the 
the role it plays in 

matic (see reviews b y iMathews fc Capriottil 
Osterbrock fc Mathewsl 119861 lOsterbrockl 



Osterbrockl 119931 and ISulentic et all |2000[ ). How 



1985 



1991 



ever, in recent years, a clear picture of t he structure a nd 
kinematics of the BLR has emerg ed (see [Gaskell 200jfor 
a detailed review) . The BLR is simply the inner exten- 
sion o f the torus dNetzer fe LaorlH993t ISuganuma et al.l 
120061: iGaskell et al.l l2007bf ) and it sha res a similar 
flatten ed distribution and covering factor (jGaskell et al.1 
I20071i hereinafter GKN). In the GKN model the BLR 
is self-shielding near the equato rial plane, s o the re is 
strong ionization stratification ([Mac Alpine 197 2]), in 
agreement with ob servat ions (see IGaskell et al.1 [2007b, 
IGaskell et all 120081 and iGaskehl 12009ft . The covering 
factor of the BLR, like that of the torus, is quite large, 
so the appearance o f the BLR is, as suggested by 
iMannucci et al.l ((199 2) . som ewhat like a " bird's nest". 
This is illustrated in Fig. 7 of lGaskehl ([2009ft . Because of 
the dusty torus, we usually see the BLR close to face-on. 
The dominant motion of the BLR is rotation, but 
there has to be significant vert ical motion to maintain 
its thickness (jOsterbrockl 11978ft . In addition to these 
motio ns there is a net inflow ([Gaskell fc Goosmannl 
2008) which is providing the feeding of the black hole. 

In sections 2 and 3 below I show that the GKN 
mode l readily explains classic "lo g arithmic" line pro- 
files ([Blumenthal fc Mathewsl 119751 : iBlddwinl 1 1975ft and 
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double-peaked, symmetr ic emission-lin e profiles expected 
from rotating disks fe.g.. lShieldslll977t ). and intermediate 
cases, without needing to invoke fundamental differences 
in the AGNs. The GKN model also reproduces observed 
transfer functions and velocity-delay diagrams. 

Despite progress in understanding the BLR, a number 
of major problems have remained. As will be discussed 
below, these include: (a) the strong asymmetry that can 
be found in the broad Balmer lines, (b) differing time 
delays of the BLR response to continuum events occur- 
ring close together in time, (c) a curious lack of corre- 
lation (or even an anti-correlation) between broad line 
variability and continuum variability over narrow veloc- 
ity ranges at certain times, (d) the contradictory and un- 
physically variable BLR kinematics implied by velocity- 
resolved reverberation mapping, (e) object-to-object dif- 
ferences in the velocity dependence of line polarization 
and, (f) changes in this polarization with time. I demon- 
strate in the remaining sections of the paper that all these 
problems and others can be readily resolved if the con- 
tinuum variability is off-axis. 

2. MODELLING BLR LINE PROFILES AND VARIABILITY 

A computer code, BL-RESP, has been developed to 
model BLR observations. BL-RESP generates many 
observable quantities and makes movies of BLR cloud 
motions seen from various viewing angles. Outputs in- 
clude profiles of emission lines, and the temporal re- 
sponse of the BLR as a function of time delay or 
"lag" (t), flux to continuum variability (the response 
is normally e xpressed as the so-called " transfer func- 
tion", *(r); iBlandford fc McKed fl98^1 . BL-RESP 



2 



GASKELL 



also gives the response as a function of projected ra- 
dial vel ocity, t>, to give so-c alled "velocity-delay dia- 
grams" (|Ulrich fc Hornelll996D . BL-RESP does no* in- 
clude detailed atomic physics. It also does not include 
the effects of scattering wh ich can modify profiles (see 
iGaskell fc Goosmannl I2008D to produ ce the blueshifting 
of high-ionization lines (Gaskcll 1982). The program can 
be used w ith the Monte Carlo radiat ive transfer code 
STOKES ((G oosmann fc Gaskelll l2007h which prod uccs 
detailed models of line and continuum polarization due 
to multiple scatterings. 

To match the distribution of gas in the GKN model, 
the BLR structure is approximated as a spherical dis- 
tribution with polar bi-cones removed. A real BLR 
probably consists of a fractal distribution of clouds 
(|Bottorff fc Ferlandll200l . but for ease of computation 
BL-RESP assumes that the BLR consists of condensa- 
tions ("clouds"). 

Because pressure is inadequate to support the observed 
thickness of the BLR there must be vertical motion of 
the gas. This is also required to explain line profiles and 
the consistency of mass estimates from the BLR. The 
inflow of material onto the black hol e is driven by the 
magn eto-rotation instability (MRI) (|Balbus fc Hawlevl 
1991). Simple energy equipartition would suggest that 
turbulence out of the plane is comparable to the MRI 
turbulence in the plane of the disc and so it is highly 
likely that magnetic instability is also responsible for the 
vertical motions. However, because there is not yet a 
satisfactory theory of these complexities, BL-RESP sim- 
ply models the vertical structure and motions by hav- 
ing clouds move in tilted circular orbits. This model 
assumes that BLR clouds are effectively collisionless and 
that their motions are dominated by gravity. This has 
the advantage of involving a known force. In reality 
magnetic fields are almost certainly playing a key role 
both in the BLR kinematics and in preventing the clouds 
from being destroyed in hypersonic collisions. Magnetic 
fields can prevent destruction of clouds in shocks be- 
cause the relevant speed in calculating the Mach number 
is the Alfven speed rath er than the sound speed (e.g., 
iKennel fc Coronitllll984l ). Such considerations are, how- 
ever, far beyond the scope of this paper. The distribution 
of tilts was taken to be uniform. 

The remaining parameters are the inner and outer 
radii, the radial distribution of clouds (taken here to be a 
power-law), the viewing angle, the distance dependence 
of the response to the continuum, and the magnitude of 
the inflow component of the cloud velocity. The models 
do not consider how local events will modify the local 
ionization structure of the BLR, although BL-RESP can 
be used to model such refinements. 

The relative sizes of the regions emitting each line and 
the effective range of radii over which that line is emitted 
are p redicted by the GKN model (see Fig. 2 of iGaskeH 
2009). The relative sizes agree w i th rev erberation map- 
ping (see Fig. 1 of IGaskell et~alll2008j ). Reverberation 
mapping can be used to set the scale factor. In this pa- 
per detailed discussion will be restricted to the Balmcr 
lines, for which the best data are available. What can 
be inferred from other lines will be discussed elsewhere. 
The relative range of radii over which the Balmer lines 
are emitted is well-determined observationally (a) from 
the widths of ^(r) recovered from reverberation map- 



ping observations and (b) from line profile fitting. These 
sizes are in excellent agreement with the predictions of 
the GKN model. For H/3, the range of radii over which 
the line is emitted is surprisingly narrow. The GKN 
model p redicts a ran ge of a factor of four or so (see 
Fig. 2 of lGaskelll[2009T ). while the best H/3 transfer func- 
tions also giv e a range of a factor of fo ur or so in radius 
(see Fig. 8 of iPiipers fc Wanders! [l99l . Such a range of 
ra dius is consistent with inner and outer radii deduced 
bv lEracleous fc Halpernl (|2003l ) from fitting disk-like line 
profiles. In contrast to this, for high-ionization lines like 
C IV (to be considered elsewhere) , the inner radius is 
very small, and the range of radius is much larger (more 
than a factor of 10). 

The typical half-opening angle of the polar cone can 
be inferred (a) from the covering factor implied by the 
observed ratio of type-1 to type-2 AGNs, and (b) from 
energy-balance arguments by considering the fraction of 
cont inuum radiation repr ocessing by the torus and BLR 
fsee IGaskell eT"aT1l2007blL The covering factor deduced 
from the energy-balance arguments is different for differ- 
ent lines , but for a line with a typical covering factor of 
~ 40% (jGaskell et al.ll2007bl ). it implies a half-opening 
ang le of 60-70°. For C I V the covering factor is higher 
(see lGaskell et al. 2007b) and the half-opening angle cor- 
respondingly smaller. 

The response of BLR emission at a given location to 
continuum variability depends on several factors includ- 
ing: the ionizing flux, Fi on , reaching the gas, the amount 
of gas at the location (which depends on the gas density 
and local filling factor), and details of photoionization 
and recombination. The ionizing flux reaching the gas 
will fall off with the distance the radiation has to travel, 
r, at least as fast as r~ 2 . It will fall off faster if ab- 
sorption and scattering along the line of sight are signif- 
icant. In BL-RESP the fall-off in Fi on is parameterized 
as Fi on oc r~ a . For the models in this paper the index, 
a, was taken to be 2. For symmetric illumination the 
effect of this distance dependence of F ion is coupled with 
the effect of radial fall off in surface density of the gas. 
This gives an overall emissivity power-law index, £. This 
was taken to be 3, t he approximate value s ugges ted by 
disk profile fitting by lEracleous fc Halpernl (|2003l ). The 
results presented here are not sensitive to the precise val- 
ues of a and £. 

In this paper I use the models only (a) to show how the 
GKN BLR model reproduces the full-range of observed 
symmetric profiles, and (b) to demonstrate qualitatively 
how off-axis energy generation solves a wide range of ma- 
jor BLR puzzles. I will not attempt to optimize param- 
eters to match specific observations. However, as will 
be shown elsewhere, it is straight forward to match ob- 
served line profiles, transfer functions, and velocity-lag 
diagrams, and in doing so to investigate constraints on 
the above parameters. 

3. BLR LINE PROFILES ARISING FROM SYMMETRIC 
CENTRAL ILLUMINATION 

The Bal mer line s of the majority of AGNs 
studied by Bald win! (|1975|) showed symmetric pro- 
files th at could be well fit by a "logarithmic" 
profile ()Blumenthal fc Mathewi I1975D . On the 
other hand, other AGNs sho w disk-like Balmer 
line profiles (see for example lEracleous fc Halpernl 
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Fig. 1. — The dot-dashed line shows the profile that would arise 
from the GKN model with parameters as described in the text if 
it were centrally illuminated and viewed from 30° off axis. The 
profile with short dashes is the same model viewed from face-on. 
The thin line superimp o sed o n this is the logarithmic profile of 
IBlumenthal fc Mathewsl lfl975h . 



1994 ; IStrateva et al.ll2003HGezari Halpern. fe Eracleousl 
2007h " lEracleous fc Halpernl (|1994l ) suggested that 



AGNs with disk-like Balmer profiles are fundamentally 
different from AGNs with strong, symmetric, centrally- 
peaked profiles. Fig. 1 shows the profiles arising from 
the GKN model with the basic parameters as described 
in the previous section. When the BLR model is viewed 
from 30° off axis, the line profile shows the expected char- 
acteristic double humps. This profile agrees well with the 
variable part of NGC 554 8 H/3 profile shown in Fig. 4 of 
IShapovalova et al.l (|2004t) . When the same BLR model is 
viewed from face-on it can be seen from Fig. 1 that a clas- 
sical "logarithmic" profile is seen instead. This demon- 
strates that disk-like line profiles can arise from the same 
sorts of BLRs that produce logarithmic profiles and thus 
supports the idea that A GNs showing disk-like and loga- 
rithmic line profiles are fundamentally the same with the 
only difference being the viewing angle. Other apparent 
differences between AGNs with disk-like and logarithmic 
profiles can also be expla i ned b y differing orientation. In 
particular, iGaskell et al.l ()2004f ) have demonstrated that 
the apparent continuum differences in AGNs seen off axis 
are simply due to reddening. 

4. OFF-AXIS CONTINUUM VARIABILITY 

The bulk of the energy in AGNs comes out in 
the t hermal "big blue bump" (see Fig. 2 of iGaskelll 
2008). This energy is produced by dissipative heat- 
ing of matter accreting o nto the black h ole. This 
was discussed in de t ail by iPring le fc Reesl (Il97l and 
IShakura fc Sunvaevl (|1973f ) in the context of optically- 
thick, but geometrically-thin, accretion disks. With only 
slight modification (the addition of heat flow from the 
inner hotter regions to the cooler outer regions) this 
scenario (or quasi-spherical accretion) reproduces the 
time -averaged spec tral energy distribution of AGNs well 
(see IGaskelll [2008). The huge problem with steady- 
accretion models, however, is that, after subtraction of 
the host galaxy starlight contribution, the amplitude of 
variability of the th ermal big blue bump is enormous 
(|GaskellH2007l l2008ft . Rather than AGN variability be- 
ing a small perturbation of the standard disk model, 
the variability mechanism is the main source of energy 



(|Gaskell fc Klimekll200l. 

It is pointed out in IGaskell (|2008j l that this highly vari- 
able energy generation must necessarily be non-axially 
symmetric. As is well known, for an accretion disk, 
the strongest contribution to the emission at a given 
frequency arises at a radius that depends on the ra- 
dial temperature profile of the disk. For a given spec- 
tral region the emission thus arises predominantly from 
within an annulus. Because information cannot propa- 
gate around a ring in an accretion disk on the timescale 
of AGN variability (which is often as short as the light- 
crossing timescale), it is clearly impossible to get an en- 
tire annulus to vary rapidly, simultaneously, and sym- 
metrically in brightness by an order of magnitude. In- 
stead, strong variability must arise in azimuthally- and 
radially-localiz ed regions. I t will thus be highly non- 
axisymmetric. iGaskell (l2008f ) argued that the complexity 



found in multi- wavelength va riability monitoring (see, for 
example, iBreedt et al.1 12009) is consistent with this pic- 
ture and proposed a number of tests of the off-axis energy 
generation scenario. 

The regions of the disk with an average temperature 
high enough to produce ionizing p hotons are su rprisingly 
close to the BLR (see Table 1 of IGaskelll l2008t ). When 
there is a flare it will be occurring in a lower-temperature 
region which will be further out and hence closer to the 
B aimer-line-emitting region of the BLR or inside it. In 
the following sections I consider the effects of this. 

5. VARIABLE LAGS 

As is well known, BLR variability lags contin- 
uum variability and t he lag gives an indicat i on of 
the size of the BLR (iLvutvi fc Cherepashchukl 119721: 
ICherepashchuk fc Lvutvilll973T j. The simplest and most 
widely used way of calculating the lag i s to cross-correlate 
line and continuum light cu r ves (IGaskell fc Sparkel 
fl98i IGaskell fc Petersonlll987h . iNetzer fc Maod (119901 ) 
pointed out that different continuum events in the inten- 
sive 1989 monitoring of NGC 5548 (|Clavel et al.1 Il99l 
iPeterson et alj 119911 ) gave lags that differed by more 
than the well-understood obs ervational er rors. This is 
discussed more extensively in iMaozl (fl99l . NGC 4151 
shows a similar discrepancy (|Malkov et alJH99 7V The 
differences in lag (a factor of ~ 50% in NGC 5548) mean 
that a single transfer function will not be able to repro- 
duce the line-flux variability corresponding to the con- 
tinuum variability. Such discrepancies are found for es- 
sentially all reverberation-m apped AGNs . 

It has been pointed out (Gaskcll 2008) that off-axis 
continuum variability readily explain differences in lags 
for different continuum events. For a BLR tilted to our 
line of sight we see the BLR responding with less delay to 
a continuum flare on the far side of the black hole than 
to a flare on the near side of the black hole. If the BLR is 
nearly face on to us, there is a similar effect if the flare is 
out of the plane of the BLR and towards us versus being 
behind the plane. 

The magnitude of these effects was modeled using BL- 
RESP. For a line-emitting region spanning a factor of 
four in radius, with a half-opening angle of 65°, tipped 
30° relative to our line of sight (conditions approximating 
the H/3-emitting region of the BLR of NGC 5548), a flare 
at the inner edge of the BLR in the equatorial plane on 
the side nearest the observer produces a lag 50% longer 
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than a similar flare on the opposite side. The difference 
increases to a factor of two for an edge-on BLR. 
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Fig. 2. — Observed transfer functions, ^(t), as a function of 
lag, r, compared with <P(t) calcu lated by BL-RESP. The fi lled 
black circles are *S?(t) derived by Piipcrs & Wanders (1994) for 
NGC 5548 for the 1990 observing season. The solid black curve 
is \P(t) for a centrally-illuminated model based on the GKN ge- 
ometry. The half-opening angle has been taken to be 65° and the 
H/3-emitting region to cover a factor of three in radius. The sys- 
tem is viewed from 20° off-axis. The long-dashed red curve shows 
VP(t) for the same BLR illuminated from the inner edge of the H/3- 
emitting region nearest the observer (at inferior conjunction). The 
short-dashed blue curve shows ^"(t) when the source of variable 
illumination is at superior conjunction, and the thin green curve 
shows ^"(t) when the illumination is at either maximum elongation. 
The op en red circles are the obse rved ^(r) for the 1989 observing 
season HPijpers fc W anders 1 991) . 

6. VARIABLE TRANSFER FUNCTIONS 

^(t) for H/? is usu ally consistent with being zero at 
zero time delay (e.g., IPiipers fc Wanders! I1994T ). As is 
expected for central illumination of the GKN model, this 
means that there is little or no gas close to the line of 
sight. The observed transfer function for NGC 5548 in 
the 1990 observing season is shown in Fig. 2. It can be 
seen that ^(r) calculated for central illumination of the 
GKN model is a good fit. Fig. 2 also shows ^(r) for the 
same BLR geometry illuminated by an off-axis variable 
continuum source, ^(t) is shown for the source being at 
various azimuths. It can be see that when the source of 
continuum variability is at inferior conjunction: (a) the 
lag is smaller (as discussed in the previous section) and 
(b) ^(t) is predicted to be non-zero at r = 0. When 
the continuum source is at superior conjunction the lag 
is longer and ^(t) is narrower. When the source is near 
either maximum elongation the difference in ^(t) from 
the central-illumination case is smaller. 

As the variable source of ionization changes position 
one expects, one therefore expects that the shape of 
^(t) will change. The only object sufficiently well ob- 
served to test this i s NGC 5548. For the years for which 
IPiipers fc~ Wandcrs (|1994[ ) give ^(t) it is similar to the 
1990 one shown in Fig. 2. The exception is the 1989 ob- 
serving season where ^(t) is non-zero at t = 1 unlike 
the other years 1 . 

1 Note that the off-axis illumination model predicts that ^(r) is 
poorly determined for 1989 because, as noted in Section 5, there 
were active regions in different locations producing different lags. 
One single ^(t) will not fit the entire observing season. As ex- 
plained in Section 13.1, in this situation programs to derive ^(t) 
attempt to improve the fit to the light curves by adding in spurious 
structure to \P(t) at large r. ^(t) is therefore uncertain at large r. 
The quoted error bars for 1989 are underestimated because they 
do not take the change in illumination position into account 



OFF- AXIS ENERGY GENERATION AND THE BROAD-LINE REGION 



7. LINE PROFILES 
Disk-like, double-peaked profiles are certainly seen 
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and reference s ther ein) . Gaskell fc Sneddenl 
iPopovic et all (|2004l ). and Bon et al.l (|2006l ) have ar- 
gued that such profiles are probably present in all BLRs 
but are simply hard to recognize when the disk is 
near to face-on. However, a big problem for simple 
disk BLR models is that usuall y one peak is signifi- 
cantly stronger than the other (IGaskelll 119831 Il996bl: 
IGezari. Halpern. fc Eracleousl |20~07|) . These line-profile 
asymmetries have been attributed to asymmetries in the 
BLR disks. Suggested causes of asymmet ries have in- 
clude d "hot spots" i n the BLR emission (|Zheng et al.l 
|1991[) , elliptical disks ( Eracleous et al.ll995T) , and warped 



disks (|Wu et al.l 120081 ). Although such models can ex- 
plain some line profiles, there are at least three ma- 
jor problems with these explanations. The first is 
that Keplerian shear will quickly wipe out asymmet- 
ric distributions of BLR gas. The second is that the 
line asymmetries can be very large. Extrem e exam- 
ples can be seen in F ig. 1 of IGaskelll (|1983[ ) and in 
IBoroson fc Lauerl (I2009D. A third problem i s that profiles 
chang e with time dVeilleux fc Zhengl Il99lt IZheng et al.l 
1991). These chan ges can be huge. For example, Mrk 668 
(= O Q 208) had a strong red p eak in the 1978 spec- 
tra of lOsterbrock fc Coheni (|1979D . the 1 982 spectrum of 
IGaskelll (|1983l ). and 1985-1991 spectra of IMarziani et al.l 
( 19931) . but it had a single-peaked Ha profile in the 1998 
spectra of IGezari. Halpern. fc Eracleousl (|2007t ) . Such 
changes are merely an extremum of BLR line-profile vari- 
ability. Line-profile changes, which can appear in periods 
as short as a week but can persist for years, are not ob- 
viously correlated with the continuum or emission-line 
flux level, and are not reverberation effects (IZheng et al.l 
119911: IWanders fc Peterson! [19971 iPeterson et al.lll999l ). 
Different parts of line profiles show different responses 
to observed continuum variability (IZheng et al.l 119911 : 
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Sergeev et al 



[2000L[200ll 12001 lShaDOvalovairaTl 2004) . 



(1999) pointed out that the changes they 



observed in the Balmer-line profiles of Mrk 6 shapes can- 
not be caused either by matter redistribution or light- 
travel-time effects. They therefore made the important 
suggestion that they are probably caused by changes in 
the anisotropy of the ionizing continuum. 

The effects on line profiles o f off-a xis illumination of 
the sort described in IGaskelll ()2008l ) can easily be in- 
vestigated with BL-RESP. Fig. 2 shows how even SDSS 
J0946+0139, the strongest and clearest case of an anoma- 
lous H/3 profile out of the ~ 9,800 l ow-red shift SDSS 
AGNs analyzed by IBoroson fc Lauerl (|2010| ). can read- 
ily be explained with off-axis illumination of the GKN 
model with the same parameters described above. The 
only difference between the model spectra in Figs. 1 and 
2 is that in Fig. 2 there is additional ionizing radiation 
being input at the inner edge of the H/3 emitting region. 
This continuum "flare", which has approximately three 
times the luminosity of the rest of the disk, is taken to 
be at maximum elongation on the receding side of the 
BLR disk. 
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Fig. 3.— The observed SDSS H/3 profile of SDSS J0946+0139 
(dots) compared with the off-axis illumination model described in 
the text (thick solid line). The NLR contribution to H/3 to the 
observed spectum has not been subtracted. The narrow emission 
line at +6000 km s" 1 is [OIII] A 4959. The dot-dashed line at 
the bottom is the same off-axis profile shown in Fig. 1 (i.e., with 
central illumination and viewed from 30° off axis). The thick solid 
line shows the effect of illuminating the GKN BLR from an active 
region located at the inner edge of the H/3-emitting region. 

8. LINE-PROFILE VARIABILITY 

As noted above, line-profile variability presents diffi- 
culties for existing models of line-profile asymmetries. 
In the off-axis variability model one might naively ex- 
pect the maximum of line-profile variability to be at the 
wavelength of the strongest displaced peak, but is usually 
observed instead to b e at the edge furthest away from the 
line center (see, e.g.. IMarziani et "all 1 19931 ). This is due 
to a combination of two effects: a change in the lumi- 
nosity of the off-axis, continuum-emitting region, and/or 
a change in position (and hence radial velocity) of the 
region. Fig. 4 shows the effects of an active region or- 
biting close to the inner edge of the Ha emitting region. 
The basic BLR parameters are as above and the motion 
of the active region is as described in the figure caption. 
From Fig. 4 one can see that adding the off-axis contin- 
uum emission shifts the edge of the blue peak bluewards 
compared to the symmetric-illumination case. The red 
peak is also shifted slightly to the blue. If the off-axis 
emission is combined with more symmetric illumination 
then the amount of shifting of the blue edge depends on 
the relative intensity of the active region and its orbital 
phase. 

Since observers often report difference profiles or 
rms profiles, Fig. 4 shows the difference profile for 
the two different epochs. Although no attempt has 
been made to match any particular observations in 
detail, Fig. 4 qualitatively matches the profile vari- 
ability seen in well-studied AGNs showing disk-like or 
highly asymm etric Balmer-line profiles. The ex ten- 
sive study by IGezari. Halpern. fc Eracleousl (|2007l ) of 
long-term profile variability in such objects provides 
many examples. For example, Fig. 4 can be com- 
pared with t he 3C 227 spectra and difference sp ectra 
in Fig. 32 of IGezari. Halpern. fc Eracleousl ([2007), and 
with the spectra and difference spectra for Mrk 668 
shown in their Fig. 30 (if the sign of the radial veloc- 
ity is reversed.) Many further examples can be seen 
in the spectra of AGNs shown in the Appendix of 
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Fig. 4. — Profile changes due to off-axis emission at the inner edge 
of the Ha-emitting region of the BLR. The thin continuous line 
at the bottom shows the profile expected for central illumination. 
The thick solid line at the top shows the profile expected when the 
BLR is illuminated close to its inner edge with the active region 
at maximum elongation on the approaching side of the BLR. The 
dashed line just below it shows how the profile changes if the active 
region moves 30° around the orbit. The dotted line at the lower left 
shows the difference between the two profiles at these two epochs. 



iGezari. Halpern. k Eracleoui (|200l . PKS 0235+023, 
for example, showed a symmetric, classic, disk-like pro- 
file in 1991, but the blue peak became much stronger by 
1998 and, as predicted here, the blue peak shifted to a 
more negative radial velocity. 

In the off-axis- variability model we expect different re- 
gions to become active at different times. If a region 
remains active long enough it will change its position 
because of orbital motion. The activity could also prop- 
agate inwards or outwards in the disk. Orbital motion 
of regions of ex c ess Ba lmer line emission was suggested 
bv IZheng etaLl (1199 1) as the cause of line-profile vari- 
ability in 3C 390.3. iGaskelll (|1996a[ ) showed that blue 
peak in 3C 390.3 showed an almost linear decrease in 
bluesh ift over a period of 20 years, but this did not con- 
tinue (jEracleous et al.l 11997). Similar apparent orbital 
motion on a shorter timescale ca n be seen in Arp 102B 
(see Fig. 7 of ISergeev et all 120001) and later in 3C 390.3 
(s ee Fig. 7 of ISergeev et al.ll2003i 

iNewman et"aT1 (|1997t l and ISergeev et all £2000) inter- 
preted the changing excesses of emission in the Balmer- 
line profiles of Arp 102B as "hot spots" orbiting with a 
period of a little over 2 years. There are several seri- 
ous problems with this. The first is that a second "hot 
spot" appearing in Arp 102B from 2000-2005 gave a 
mass that differed by ~ 30 times t he measuring error 
(jGezari. Halpern. k Eracleousl 12007) . A second prob- 
lem is that for other AGNs with disk-like profiles the 
timescales are totally inconsistent with the black hole 
masses estimated from the relationship b etween stel- 
lar ve locity dispersion and black hole mass ()Lewis et al.l 
2010). A third problem is that while excess flux in a 
line profile is observed at times to change between the 
red and blue sides of a line (as is expected if it is due to 
a region of excess line emission orbiting), the observed 
radial velocity curve can be incompatible with simple or- 
bital motion. Fig. 5 shows an example of this. 

In the off-axis variability model the excess emission is 
due to temporarily active regions preferentially illumi- 
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Fig. 5. — Radi al velocity curves for regions of ex cess Ha emission 
identified bv lGezari. Halpern. feE raclcous (2007) in Mrk 668. The 
horizontal lines are the average radial velocities for the periods 
indicated. 

nating BLR gas at a particular velocity. An apparent 
change in radial velocity of excess emission, such as the 
changes shown in Fig. 5, is due to a region at one location 
becoming inactive and another region becoming active in 
a different location. For high-mass black hole systems the 
orbital period will be long, so we expect that the radial 
velocity will not change slowly as an active region moves, 
but will change suddenly. This sudden change in velocity 
seems to be what is happening in Mrk 668 (see Fig. 5). 
In this case what matters is how the active regions pro- 
ducing the ionizing continuum radiation change, not how 
the structure of the BLR gas is changing. 

9. VELOCITY-DEPENDENT FLUX CORRELATIONS 

In their studies of profil e variability o f Arp 102B , 
NGC 4151, and 3C 390.3 ISergeev et all (l2000l l200l 
2002) found three curious things in all three objects that 
have hitherto not been satisfactorily explained: 

1. Ha showed strong rms variations in different ve- 
locity bins while the total line flux changed only a 
little. 

2. The integrated line flux has a much better correla- 
tion with continuum flux than any individual ve- 
locity bin does. 

3. Relatively narrow velocity regions of line profiles 
can at times show a strikingly weaker correlation 
with the continuum (or even anti-correlations.) 

All three of these effects can readily be explained by the 
off-axis- variability model. Of the three the most remark- 
able is the lack of correlation or even anti-correlation 
with the continuum over re latively small veloci ty inter- 
vals (see especially Fig. 9 of ISergeev et alj|2001l ) Fig. 6 
illustrates where a narrow range of radial velocity has to 
come from. For illustration purposes the BLR disk has 
been taken to be infinitely thin - for a more typical disk 
scale height the region where a narrow range of velocity 
is coming from is more diffuse. Different velocities in a 
Balmer line profile correspond to different distances of 
the centers of the crescents along the projected major 
axis of the BLR. The narrower the velocity region, the 
narrower the crescent. 

The first of the three effects noted by Sergeev et al. 
is readily explained: many different off-axis regions are 
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Fig. 6. — Schematic illustration of the Balmer-line-emitting re- 
gion of an idealized, geometrically-thin BLR disk viewed from 30° 
off-axis. The crescent-shaped region on the left shows the gas emit- 
ting in an arbitrary narrow range of observed radial velocity. In a 
realistic thick BLR the boundaries of the region are not sharp. 



contributing to the variability; each individual region will 
predominantly influence the line only over a narrow range 
of radial velocity; but the integrated flux of the line is 
the sum of the contributions of all the regions. The lower 
rms variability of the line as a whole compared with indi- 
vidual velocity bins is thus simply the result of averaging 
over more events. 

The explanation of the second effect is similar: the 
observed continuum variability is the sum of all the indi- 
vidual off-axis events and is thus better correlated with 
the total line flux (which also depends on the sum of all 
the individual events) than with the line flux in an indi- 
vidual velocity bin which is dominated by the variability 
of a subset of the events. 

The third effect arises because a minor off-axis event 
far out in the BLR disk might not make a detectable 
contribution to the total continuum flux variability seen 
by the observer, but it will have a strong effect on the 
localized region it is close to in the BLR. At that loca- 
tion the variability of the closest off-axis-variability re- 
gion can easily dominate over the combined variability of 
more distant regions. If the local variability happens to 
be out of phase with the variability of the brighter, more 
distant, continuum-emitting regions which dominate the 
observed continuum light curve, the result will be to give 
an anti-correlation between the continuum flux and the 
line flux in the velocity range the localized region is con- 
tributing to. Conversely, if the variations are in phase 
the correlation will be enhanced over the velocity range. 
Of these two effects, weakening the correlation will be the 
more obvious effect if the correlation coefficient is plot- 
ted against velocity. Strengthening the correlation over 
a velocity range will be more obvious if the significance 
level is plotted against velocity (e.g., as in iLewis et al.l 
120101) . 
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Fig. 7. — The effect of off-axis variability on the velocity- 
dependent correlation coefficients. Continuum variability has sim- 
ply been assumed to be the result of events of different amplitudes 
at different radii. Approximately half a dozen regions were taken 
to be active at any moment and the sampling time of the observa- 
tions was taken to be longer than the lifetime of the events. The 
continuum flux at any instant is the sum of fluxes of all the indi- 
vidual events. Each event was assumed to most strongly influence 
a velocity range of ±1000 km s _1 . The distribution of events was 
taken to be uniform in line-of-sight velocity. Different lines corre- 
spond to different seeds of the random number generator and show 
the range of curves seen. The correlation coefficient between the 
continuum and the line as a whole was always greater than 0.9. 

The effect of off-axis variability on the velocity- 
dependent correlation coefficients has been qualitatively 
modeled in Fig. 7 (see figure caption for details). Note 
that this is to illustrate how strong variations in the cor- 
relation coefficient can readily be produced and that no 
attempt has been made to input a realistic distribution 
of continuum events or to match the observations of any 
object. 

The main thing to note from Fig 7 is that the off- 
axis-variability model readily produces narrow velocity 
regions where the line flux is not correlate d with the con- 
tinuum . If we compare Fig 7 with Fig. 9 of lSergeev et al.l 
we can see that in NGC 4151 the velocity ranges 
of non-correlation with the continuum flux can be even 
narrower than produced by the off-axis- variability model 
in Fig. 7 (in 1997-1998, for example). Similarly narrow 
regions of enh anced correlation can be found in CBS 74 
(see Fig. 1 2 oflLewis et al.lj201Ct ) and PKS 0921-123 (see 
Fig. 13 of lLewis et alJl2~010( ). The narrowness of these 
regions implies that the region of the BLR influenced by 
an off-axis event can be quite small. 

Three importa nt things that analy zes such as those 
shown in Fig. 9 of lSergeev et al.l <|2001f > do is to establish 
that: 

1 . there is significant off-axis ionizing continuum vari- 
ability occurring out at the radius of the Balmer- 
line-emitting region, and 

2. the off-axis variability is not far above or below the 
BLR disk. 

3. the Balmer line emitting region is indeed in a flat 
disk. 

This is because (a) the radii are accurately pinned down 
by the velocities at which poor line-continuum correla- 
tions are seen, and (b) neither a very thick BLR nor con- 
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tinuum events far from the equatorial plane will produce 
the sharp velocity dependence of the correlation coeffi- 
cient. 

Note that, in general, the spectral region correlated 
best with the continuum is not the region that shows 
the strongest variability, although in some cases it can 
be (e.g., at —750 km s _1 in Pictor A - see Fig. 11 of 
iLewis et al.ll2010l) . What matters is the phase of the local 
variability. For the most variable part of a line this might 
or might not be in phase with the dominant continuum 
variability. 

10. VELOCITY-RESOLVED REVERBERATION MAPPING 

Velocity-resolved reverberation mapping shows that 
the motions of BLR gas clouds are gravita tionally domi- 
nated rather than arising from an outflow (|Gaskelllll988l ; 
iKoratkar fc Gaskellll989D . This is because the red and 
blue wings vary closely in phase on a light-crossing 
timescale. On average the red wing leads slightly, 
imply ing a slight net inflow (see Gask ell fc Goosmannl 
2008 for a summary of the case for an inflow com- 
ponent of motion). While the majority of velocity- 
resolved reverberation mapping observations a re con- 
sistent with this picture (see iGaskelll l2010bD . there 
are cases where velocity-resolved reverberation map- 
ping clearly shows a strong signature of apparent out - 
flow (IKollatschnv fc DietricIj|1996l : lDennev et al]|2009bl ). 
IDennev et alj (|2009br i suggest that differing kinematic 
signatures reflect real object-to-o bject di f ference s (e.g., 
a strong outflow in some AGNs). IGaskelll (l2010bTl . how- 
ever, points out that differing kinematic signatures have 
been seen in the same objec t s at different time s (e.g., 
IKollatschnv fc Dietrich! fl99l iWelsh et all 120071 ') The 
most spectacular exa mple occurred in t he 1989 moni- 
torin g of N GC 5548 (IClavel et al.lll99lt iPeterson et al.l 
119911) where IKollatschnv fc Dietrich! (|1996f) found strong 
outflow and inflow signatures for C IV and the B aimer 
lines in successive outbursts. Since these comparable 
outbursts were separated by only 100 days (about the 
light-crossing time for the outer BLR) a real change of 
direction of motion of the entire BLR is out of the ques- 
tion. 

Off-axis illumination provides a simple explanation of 
what is going on. If from our viewpoint the off-axis con- 
tinuum variability is displaced towards the approaching 
side of BLR, the blueshifted side of the line profile will 
vary before the redshifted side. This gives an apparent 
"outflow" signature. Conversely, if the continuum vari- 
ability is displaced towards the receding side, the red 
wing variability leads the blue to give an "inflow" signa- 
ture. The velocity dependence of the lags is similar in 
both cases; only the sign of the velocity is different. The 
solid curve in Fig. 8 shows the velocity dependence of the 
H/3 lag predicted by BL-RESP for a 4:1 ratio of inner to 
outer radius, a 65° half opening angle, a 30° tilt to our 
line of sight, and a "flare" at maximum elongation right 
at the inner edge of the BLR. The velocity dependence 
of the lags was also calculated for an accelerating wind 
and is shown as a dotted line in Fig. 7. These theoretical 
curves can be compared with the four best examples of 
a strong velocity depende nce in the K(3 lag: the 1993- 
1995 monitoring of Mrk 6 (|Sergeev et alJll999ri. the 2007 
monito ring of NGC 3516 and NGC 3227 (IDennev et al ' 
l2009bl) . and the 2008 monitoring of Mrk 40 ()Bentz et~a! 
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Fig. 8. — Observed radial velocity dependence of lags in Mrk 6 
(solid circles), NGC 3516 (squares) , NGC 3227 (open circles joined 
by dotted line), and Mrk 40 (triangles). The velocity scale for 
NGC 3227 has been reversed. The solid line shows the predicted 
lags as a function of radial velocity for the GKN model with off- 
axis continuum variability. The dashed line shows the predicted 
lags for an accelerating outflowing wind. 



I2008f) . NGC 3227 shows the opposite sign of the velocity 
dependence to the other three, so it has been plotted with 
the velocity scale reversed. The four AGNs have different 
line widths and total lags. To facilitate comparisons the 
lags of both observations and models have been normal- 
ized to a 15-day lag at zero velocity. Choosing an appro- 
priate normalization for the line widths is harder since 
different commonly used line width estimators give sig- 
nificantly different widths. The velocity normalizations 
have therefore been adjusted to minimize residuals from 
the off-axis flare model (the thick solid line). Error bars 
in the lags have been omitted for clarity, but they can 
be found in the original papers. They are comparable to 
the scatter about the theoretical line at low velocities, 
and larger at the highest velocities. 

The first thing to notice in Fig. 8 is that an accelerat- 
ing outflowing wind model is a poor fit to the NGC 3227 
lags. In the wind model the high-velocity red wing should 
show the largest lag, but the observed lag decreases in 
the wings. It can be seen instead that, after reversing the 
velocities, the velocity dependence of the NGC 3227 lags 
is completely consistent with the velocity dependence for 
the other three AGNs. This supports the idea of a com- 
mon origin of the velocity dependence of the lags. Fi- 
nally, it can be seen that there is good agreement with the 
velocity dependence predicted by the off-axis-variability 
model. Other than the scaling in velocity, there has been 
no attempt to optimize model parameters. The shape 
of the theoretical curve is largely independent of these. 
The main effect of changing parameters is to increase or 
decrease the amplitude of the blue/red difference. For 
example, a greater tilt to the line of sight increases the 
difference, while moving the "flare" closer to the black 
hole, or changing to phase away from maximum elonga- 
tion reduces the difference. The off-axis variability shown 
in Fig. 8 produces about the maximum effect possible for 
the viewing angles expected for most AGNs. This should 
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not cause concern because the four AGNs were chosen to 
be extreme cases. 

Although an outflow (either accelerating or constant- 
velocity) is excluded even for NGC 3227, infall is not 
excluded for the other three. To reproduce the strong 
red-wing/blue-wing diffe rence the inflow w ould need to 
be close to radial freefall (|Bentz et al.ll2010l have indepen- 
dently obtained the same result for Mrk 40 = Arp 151). 
BL-RESP shows that the velocity dependence of the lag 
for pure infall is essentially indistinguishable from the 
off-axis model shown in Fig. 8. However, pure radial 
infall does not explain how different kinematic signa- 
tures are seen at different times. Mrk 6, for example, 
showed a much les s marked red-wing/ blue-wing differ- 
ence m 1996-1997 (|Sergeev et al.lll999j ). More extreme 
examples of changing kinematic signatures have been 
mentioned above. 

11. HYDROGEN AND HELIUM LINE RATIOS 

Although local BLR ionization perturbations caused 
by nearby continuum variability have not been included 
here there will be such effects. In the GKN model 
high-ionization emission (e.g., Hell) comes from much 
closer in than the B aimer line emission and hence is seen 
at higher velocities. Off-axis variability, if it produces 
enough flux with energies above 4 Rydbergs, will tem- 
porarily cause additional enhanced He II emission further 
out in the disk (and hence at lower velocities) . The model 
therefore predicts that Hell, like the Balmer lines will 
show anomalous variability behavior over the same nar- 
row velocity ranges and at the same time as the B almer 
lines. This is precisely what ISergeev et al.1 (|2001[ ) find 
(see their Fig. 9). 

It has l ong b een known empirically (jShuderl 119821 : 
I Crenshaw! 119861 ) that the Balmer decrement is flatter 
for higher velocity BLR gas. Photoionization modelling 
shows that the Balmer decrement flattens (H/3/Ha in- 
creases) with increasing io nizing photon flux (|Kwanlll984t 
iSnedden fc Gaskelll l2007h . In the off- axis illumination 
model the photon flux is higher in the regions closest to 
the continuum sources. We therefore expect that these 
regions will have flatter Balmer decrements. This means 
that "bumps" in the line profiles due to off-axis illumi- 
nation will be more prominent in H/3 than in Ha. In- 
spection of Fig. 2 of IShapovalova et al.l (|2010t ) reveals a 
number of epochs when this is the case (e.g., in July 1996, 
January 1999, and January/February 2000). 

12. LINE POLARIZATION 

Polariza tion gives powerf ul insights into the structure 
of AGNs (|Antonuccil l2002h . Polarization is extremely 
sensitive to departure from symmetry about the line of 
sight and therefore provides a powerful test of the off-axis 
illumination model being presented here. 

Li ne emission from a sym metrical disk will be polar- 
ized (|Chen fe Hafpern|ll990f). Spectropolarim etric obser- 
vations of Arp 102B (|Antonucci et al.lll996[) did indeed 
show polarization of t he line but this was in consistent 
with the predictions of lChen fc Halpernl Jl990). The line 
polarization was twice as great as predicted and, contrary 
to the model predictions, the polarized flux was higher 
in part of the blue peak than in the red peak, and it also 
dropped too low in the far blue wing. The polarization of 
the blue side changed dramatically with wavelength, con- 



trary to the prediction. The polarization angle changed 
relative to that of the continuum by ~ 30° compared 
with the continuum rather than the predicted ~ 90°. 

Spectropolarimetry of AGNs in general s hows a diver- 
sity o f line polarization (see examples in the ISmith et al.l 
2002 spectropolarimetric atlas). It is normal for the po- 
larization angle to vary across a line and also for a line 
to be less polarized on one side and more polarized on 
the other. Differing line polarization behaviors have been 
interpreted as consequences of different kinematics and 
scattering geometries in different objects with, for exam- 
ple, scattering from a rotating region domin ating in one 
objec t and scattering from a wind in another (jAxon et al.l 
l2008h . 

The off-axis continuum emission discussed here pro- 
duces both line and continuum polarization. The polar- 
ization behavior is easy to understand qualitatively. Po- 
larization arises from a departure from symmetry about 
the line of sight. A face-on disk, for example, has no 
asymmetry about the line of sight, and hence the inte- 
grated light has no polarization, but the polarization in- 
crea ses as the disk inclinat i on is increased (see, for exam- 
ple, [GoosmanE^G^Mi HQ03). Off- axis illumination 
similarly produces strong polarization, but the polariza- 
tion measured by the observer is diluted by the combined 
effects of illumination from other directions. As discussed 
above, the gas nearest a source of off-axis illumination is 
brightest and shows the strongest variability. However, 
because polarization measures asymmetry, the gas near- 
est an off-axis source of illumination will show the lowest 
polarization because the illumination is most symmetric 
nearest the source of illumination. BLR gas on the oppo- 
site side of the disk will show the highest polarization as 
a result of that particular off-axis continuum source even 
though the illumination on the opposite side will not be 
particularly bright. This means that if the source of off- 
axis illumination is on the receding side of the BLR, the 
red side will show the lowest polarization (because the 
illumination is relatively symmetric) , while the blue side 
of the line will show higher polarization. If one looks at 
the polarized flux, the line will appear to be blueshifted. 
If the source of off-axis illumination is on the other side 
of the line the effects are simply reversed in velocity. 

Detailed quantitative modeling needs a code such as 
STOKES which handles multiple-scatterings and the 
atomic physics of line polarization. Such modelling will 
be discussed elsewhere, but since the primary effects are 
geometric the qualitative behavior of the polarization can 
be derived from BL-RESP output. Fig. 9 shows the total 
line flux, the relative percentage of line polarization, and 
the effect on the polarization position angle of the line 
plus an arbitrary polarized continuum. For illustration 
purposes the BLR is assumed to be excited by an active 
region at an azimuth 45° from the line of nodes in the 
middle of the Balmer line emitting disk. 

The off-axis illumination model predicts that, in ad- 
dition to this characteristic velocity dependence of the 
degree of polarization, there will always be a change in 
the angle of polarization as one goes across the line cen- 
ter, except in the unlikely cases of the off-axis illumina- 
tion source lying exactly on the line of nodes or going 
through inferior or superior conjunction. In all models 
where the angle of BLR polarization is at least slightly 
different from the angle of continuum polarization there 
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Fig. 9. — Total line flux (solid black line), line polarization 
(long-dashed blue line), and combined polarization of the line plus 
continuum (short-dashed burnt orange line) for a BLR illuminated 
by off-axis emission in the middle of the BLR at 45° from the line 
of nodes. The light dotted line (burnt orange) is the continuum 
polarization position angle. For the total flux and line polarization, 
zero is at the bottom of the graph but the vertical scale factor is 
arbitrary. 

will be a change in the position angle of the polarization 
at an emission line. Because the change in angle is sim- 
ply the result of adding the Stokes vectors the change 
will be greatest at the line center. This can be seen in 
Fig. 9. 

Spectropolarimetric observations strongly support 
both the predicted degree of polarization and position- 
angle behaviors. The wavelength at which the excess flux 
is centered is indeed also the velocity at which the per- 
centage polarization is lowest. If we go to minus that 
velocity on the other side of the line the polarization is 
high. T his can be seen in Figs. 6, 7, 10, 13, 14, 19, 
and 23 of ISmith et al.l (l2002f) and it is especially clear in 
E1821+643 ([Robinson et al.ll2010j ). These examples in- 
clude four cases of the flux excess and polarization mini- 
mum being on the red side of the line, and three cases of 
them being on the blue side of the line. The velocity de- 
pendence of the degree of polarization shown in Fig. 9 of 
this paper can be directly compared with th e continuum- 
subtra cted polarization shown in Fig. 23 of ISmith et al.l 
(2002). The S-shaped change in position angle shown 
in Fig. 9 as one goes across th e line profile can b e seen 
in the majority of AGNs in the ISmith et al.l (|2002f) atlas 
(see Figs, lab, 2, 5b, 6, 7ab, 10, 11, 12ab, 13, 14, 15, 
16, 20, 22, and 23). More detailed analysis of individual 
objects will be presented elsewhere. 

The alternative explanations offered in the literature 
for the diverse spectropolarimetric observational results 
are that the wavelength dependence is the result of ro- 
tation in some objects, and scattering off a wind in oth- 
ers. (jRobinson et alJl2010t ) have also made the sugges- 
tion that in E1821+643 the the broad-line region has two 
components moving with different bulk velocities away 
from the observer and toward a scattering region at rest 
in the host galaxy because of gravitational recoil due to 
anisotropic emission of gravitational waves following the 
merger of a supermassive black hole (SMBH) binary. In 



contrast to these diverse explanations in the literature, in 
the model presented here, all of the spectropolarimetric 
observations are simply explained by one process: off- 
axis illumination. 

13. DISCUSSION 
13.1. The Limits of Reverberation Mapping 

Off-axis variability severely limits what can be 
achieved with reverberation mapping. It means that 
goals such as trying to map p ossible spiral density 
waves m BLR s (iHorne et al.ll2004D are unattainable. As 
pointed out in lGaskelll (|2010blT making further progress 
in reverberation mapping is not simply a matter of get- 
ting better sampling and higher signal-to-noise ratios; it 
is the nature of continuum variability in the AGNs and 
the response to it which is the fundamental limit. No 
matter how good the sampling, there is a strong risk of 
getting very misleading results from observing campaigns 
of short duration because the results are dominated by a 
small number of events. It is more important to have 
longer monitoring campaigns rather than denser sam- 
pling. In the past a lag has been considered to be well- 
determined if one complete event (a clear rise and fall or 
vice versa) is well sampled in both the continuum and line 
light curves. However, we have seen here that because 
of off-axis variability, lags from single events are untrust- 
worthy. Even two or three successive events could give 
a misleading result because, as we have seen, active re- 
gions can persist for much longer than the light-crossing 
time. Unfortunately the need to observe more than a sin- 
gle complete event adds to what is already a very labor 
intensive process. 

The off- axis- variability model means that we have to 
look in a new way at variability data in general and rever- 
beration mapping in particular. Events have to be looked 
at individually. Attempts to reconstruct Vl/(r) have of- 
ten looked at as long a series of observations as possible. 
Since the lines respond differently to different events, a 
single transfer function will fit some events poorly. A fit- 
ting routine that minimizes residuals by adjusting ^(r) 
will attempt to improve a poor fit to the line response 
by adding in a response to a much earlier event (i.e., at 
larger r.) For methods that force ^(r) to be positive 
(m ethods such as th e maximum entropy method used 
bv lKrolik et~aT1ll991| ) this will give spurious features at 
large r and an erroneous impression that the BLR is 
larger than it really is. Examples of such features can 
be seen in the transfer function s of CIV and Mg ll for 
NGC 5548 shown in Fig. 11 of iKrolik et all (|l99ll ). On 
close inspection it can be seen that even with the addi- 
tion of such features the line light curves systematically 
fail to follow the convolution of the observed continuum 
curve with *( T ) (e.g., around MJD 7540). 

To learn from reverberation mapping it is important 
to look at line profiles as well as integrated fluxes. For 
planning future reverberation mapping campaigns, good 
spectral resolution and high signal-to-noise ratios in lines 
are needed. A lthough the sort of "high-fidelity" mapping 
envisioned by lHorne et al.l (|2004f ) is unattainable, rever- 
beration mapping, if properly interpreted, still tells us a 
lot about the structure of the BLR and about the nature 
of continuum variability. A good illustra tion of what can 
be acc omplished is the recent study of I Jovanovic et all 
(2010) who track the movement of the dominant region 
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of off-axis variability in 3C 390.3 observed in 1995-1999. 
In agreement with what is found here, they find that the 
off-axis emission dominates over axisymmetric emission. 

13.2. Black Hole Masses 

Off-axis variability has multiple implications f or the 
estimation of black hole masses from the BLR (jDibail 
fl97l . To get masses from the BLR we need to know 
that the BLR is gravitationally bound. One then needs 
to know an effective radius, an effective velocity, and the 
appropriate scaling factor. The gravitational domina- 
tion of the BLR was e stablished wit h velocity-resolved 
reverberation mapping (Gaskcll 1988), but, as discussed 
above, some velocity-resolved reverberation map ping 
(|Kollatschnv fc Dietrich! 119961; iDennev et al.ll2009bH has 
given contradictory results. The good news is that ap- 
parent outflow signatures from velocity-resolved rever- 
beration mapping are the result of off-axis variability 
rather than outflows. The bad news is that because of 
the dominance of off-axis variability, reverberation times 
(our most direct way of getting effective radii) cannot be 
tr usted when o nly a single event is observed (see above) . 

iDibail (|1977t ) introduced a simpler method of estimat- 
ing the effective radius of the BLR from the AGN lu- 
minosity. Because of its ease of use, the Dibai method 
has now been applied to t ens of thousands of AGNs. 
iBochkarev fc Gaskelll (|2009l ) show that even the earliest 
mass estimates by the Dibai method agree well with re- 
verberation mapping and argue that this is because of 
the similarity of continuum shapes, BLR geometries, and 
BLR kinematics. This results in the tight relationship 
between BLR radius and continuum luminosity. Much 
of the scatter in this relationship must be consequence of 
off-axis variability. As the quality of reverberation map- 
ping data improve, the tightness of the relationship is 
there fore expected to improve (see Fig. 6 of lDennev et al.l 
2010). However, in the absence of adequate reverbera- 
tion mapping coverage of a given AGN it is preferable to 
estimate the black hole mass by the Dibai method. 

It has been common in recent years to estimate the 
effective velocity of BLR gas from rms spectra. The jus- 
tification has been that the width of the rms spectrum is 
a better indicator of the velocity of the variable compo- 
nent of the BLR gas than the width of average spectra 
is. The off-axis-variability model says instead that the 
rms spectrum will have lumps and spikes in it because 
of the multiple off-axis variability events and hence that 
it will be a poorer indicator of the effective velocity of 
the BLR than the average spectrum. Lumps and spikes 
are clearly observed in many rms spectra (see for ex- 
ample, the stron g blue spike in the 2007 rms spectrum 
of NGC 5548 in lDennev et al.ll201Q[ ). The main use of 
rms spectra is in removing non-varying narrow-line re- 
gion contributions. 

Assuming that the GKN model is an approximately 
correct description of the BLR, the potentially most se- 
rious problem in determining masses is the orientation 
effect. It can be seen from Fig. 1 that for the range of 
inclinations over which the BLR can be readily observed 
in the optical (0 to 30°) the observed FWHM (the most 
commonly used indicator of the effective BLR velocity) 
varies by about a factor of three for the same BLR. Since 
virial mass estimates depend on the square of the veloc- 
ity, the black hole mass is being overestimated by an or- 



der of magnitude for highly-inclined systems compared 
with face-on ones. The inclination effect thus needs to 
be allowed for in estimating masses. The magnitude of 
the effect depends on the thinness of the BLR. The best 
indicator of the inclination of the BLR is the flatness of 
the peak of the line profile after removing the narrow-line 
contribution. 

13.3. Looking for Signatures of Sub-parsec Supermassive 
Binary Black Holes 

It is now agreed that massive galaxies (ellipticals and 
disk galaxies with classical bulges) grow through merg- 
ers. The central black holes of merging galaxies will 
form binary b l ack h oles ()Begelman. Blandfo rd. fe Reei 
2000). Gaskcll (|1983ft suggested that close supermassive 
binary black holes could be detected through their ef- 
fect on the BLR, and, in particular, that displaced broad 
line peaks could be a consequence of the orbital mo- 
tions of close binary black holes. Reverberation mapping 
of 3C 390.3, however, showed that on a light-crossing 
timesc ale the red and blue peaks varied near sim ultane- 
ously ijDietrich et al.l [l998t lO'Brien et al.lll998t ). This 
strongly rules out bi nary black holes as the cause of 
the displaced peaks (Gaskell & Sncddcn 1999; Gaskcll 
l2010al) . Nonetheless, binary black holes certainly must 
form, and because of their importance in the evolu- 
tion of black holes and galaxies, increasing attention has 
been devoted to finding examples of clo se (sub-parsec 
scale) binaries. iBoroson fc Lauerl ([2009) have recently 
suggested that the asymmetric displaced BLR peaks in 
SDSS J1536+0441 are due t o a binary black hole sy stem 
(i.e., as in iGaskelll Il98l ) IChornock etail (|2010f ) and 
lOaskelll 42010aD argue that SDSS J1536+0441 is sim- 
ply an extreme example of an AGN with asymmetric, 
disk-like emission instead. The demonstration here (see 
Fig. 2) of how easily off-axis emission can produce ex- 
treme profiles (J1536+0441 is less extreme than the ex- 
ample in Fig. 2) adds strong support to J1536+0441 hav- 
ing a normal BLR. 

The recent debate over J1536+0441 does raise the in- 
teresting and important question of whether sub-parsec 
supermassive black hole binaries can be detected through 
anomalous BLR profiles. Because off-axis emission can 
readily produce complex profiles, the answer is unfortu- 
nately almost certainly "no" . 

13.4. Other Possible Causes of N on- axisymmetric 
Emission 

Although accretion disk variability must necessarily be 
off-axis and although, as discussed above, such variability 
successfully explains a wide range of hitherto puzzling 
BLR observations, there are a number of other plausible 
causes of non-axisymmctric variability. 

1. iGaskell fc Klimekl (|2003ft point out that the short 
timescale of the strong variability of AGNs requires 
the transmission of large amounts of energy in 
particles with relativistic or near-relativistic ener- 
gies. Such moti ons will result i n anisotropic emis- 
sion (beaming). IGaskell (]2006f ) proposed that the 
combination of correlated and uncorrelated short- 
timescale, multi-wavelength variability in AGNs 
could be a consequence of s uch anisotropi c high- 
energy emission (see Fig. 5 of lGaskelll 12000 .) 
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2. Anisotropy and variability can also be due to block- 
ing o f radiation fro m the inner regions of the 
AGN (|Gaskellll2010bh . There is abundant evidence 
now for variable X-ray absorption, often with par- 
tial coverage, in the inner regions of AGNs (see 
iTurner fc Milled 120091 IGaskelll l2010bl and refer- 
ences therein) . iGaskelll ( 201 Obi ) points out that this 
would attenuate radiation hitting the BLR in cer- 
tain directions. 

3. Another thing that could block radiation from 
the center is war ping of a thin accretion disk 
(jKinnev et al.ll200l . Such warping is well-known 
in Galactic X-ray binaries such a s SS 433 and 
Her X-l (see lOgilvie fc Dubusll200ll for discussion) 
and there are a number of possible mechanisms 
which could produc e a similar warping in AGNs 
()Kinnev et al.ll2001l) . iBentz et all (|2010[ ) have re- 
cently shown that an alternative explanation of the 
apparent strong inflow signature in Mrk 40 could be 
anisotropic illumination of the BLR due to warping 
of the inner accretion disk. 

13.5. Unanswered Questions 

Although the combination of the geometry of the GKN 
BLR model and off-axis variability has strong support 
and is successful in explaining many observations of line 
and continuum variability, there are many unanswered 
questions. These include: 

1. What is the mechanism producing the strong vari- 
ability? 

2. What is the radial and vertical distribution of the 
regions of strong variability ("flares")? 

3. What is the energy distribution of flares? (and how 
does it depend on radius?) 

4. What governs how long a continuum-emitting re- 
gion is active? 

5. What maintains the vertical motions of the BLR 
gas? 

6. What happens when rapidly moving regions of 
higher density BLR gas ("clouds") collide? Why 
are the clouds not destroyed? 

7. How much damage is done to the BLR disk by the 
strong variability close to it? 

Although a lot of further observational and theoreti- 
cal study is needed, there are indications of some possi- 
ble answers. For example, we can infer that the flares 
are not far above the BLR because it would be hard to 
produce the observed frequent lack of a line/continuum 
variability correlation over a narrow velocity range (see 
Section 9). If the flares are high above the BLR their 
effect would be spread out over too large a range of ve- 
locity. We can also say that although AGN variability 
is not azimuthally symmetric on short timescales, if we 
average over a long enough time it is expected to be sym- 
metric. Just how long is "long enough" is not clear. For 
some AGNs, asymmetric emission has persisted for as 
long as the AGN has been observed. 



For the radial distribution of activity, there certainly 
is some randomness, but if energy generation is av- 
eraged over a sufficiently long time period, the ra- 
dial dependence of the energy production mus t match 
that expected for a steady accretion flow (see IGaskelll 
2008). The frequency of structure in line profiles as 
a function of velocity and time can reveal information 
about the radial d is tribut ion of continuum variability. 
iFlohic fc EracleouH (|2008f ) have already deduced that 
there is an apparent increase in small-scale structure in 
the outer part of the H/3 emitting region. They inter- 
pret this as increased clumping in the outer BLR due 
to gravitational instabilities. In the off-axis variability 
scenario this structure is due to small scale variability in 
the continuum rather than clumping of the BLR. 

13.6. Predictions 

The off-axis illumination model presented here makes 
many testable predictions. The strongest of these, and 
the easiest to test, are the polarization predictions since 
polarization is highly sensitive to departures from sym- 
metry about the line of sight. 

As discussed in Section 12, the off- axis illumination 
model predicts that there will be a minimum in the po- 
larization near the location of the continuum source. The 
velocity of excess emission in a line profile will be the 
velocity of the minimum in the line polarization. This 
can be checked by looking at the polarization of AGNs 
showing asymmetric broad lines. In contrast to the pre- 
dictions of the off-axis illumination model, if a warp in 
the disk causes only part of t he outer disk to be centrally 
illuminated (as suggested bv lBentz et al.ll2010l ). the illu- 
minated side of the disk will be highly polarized and the 
polarization will peak at the peak of the line emission. 
This would give a quite different wavelength dependence 
of the polarization from that predicted by the off-axis 
illumination model. 

The most important tests come from polarization vari- 
ability. As different regions in different locations become 
active the general level of polarization and its velocity 
dependence will change. Changes in the velocity depen- 
dence of polari zation have a lready been observed (e.g., 
in NGC 4151; iMartell [l998h ■ The off-axis-illumination 
model predicts that the velocity dependence of the per- 
centage of polarization will reverse when the excess flux 
in a line changes sides. If there is a small change in the 
position of the main region of illumination due to orbital 
motion or other drift of the active region, this will show 
up as a change in position angle of the polarization and 
in the velocity dependence of this angle. The best can- 
didates for testing these prediction are AGNs where the 
orbital timescale is not too long and there has been a 
history of the relative intensities of red and blue excesses 
switching. 

Light echoes in broad -band polarized flux have al- 
ready been reported (see lGaskell et aLll2007a|) . Measur- 
ing velocity-dependent spectropolarimetric lags should 
be a powerful tool for distinguishing between models of 
the structure and kinematics of the BLR and continuum. 
Although it is beyond the scope of the present paper, it is 
straight forward to use BL-RESP and STOKES to pre- 
dict time-dependent and velocity-dependent spectropo- 
larimetric variations for the off-axis variability model and 
other models. 



OFF- AXIS ENERGY GENERATION AND THE BROAD-LINE REGION 



13 



The off-axis variability model makes predictions about 
the velocity dependence and variability of line ratios (see 
Section 10). Bumps due to off-axis illumination should 
be stronger in H/3 than in Ha. He I and Hell should 
show the same anomalous responses over the same nar- 
row velocity ranges as the hydrogen lines. 

In checking the model presented here there is obviously 
a need to make everything consistent (e.g., changing lags, 
changing velocity dependence, line profiles, spectropo- 
larimetry). Note though that things will not always be 
quite as simple as might naively be expected. For exam- 
ple, the brightest region of off-axis-illumination might 
not be the most variable during a given time period. 

14. CONCLUSIONS 

Although we do not understand the physical mecha- 
nism causing AGN variability, it must be highly non- 
axisymmetric. This is a major paradigm shift in our 
under standing of AGNs. It has been argued in iGaskelll 
(2008) that off-axis continuum emission explains many 
puzzles of AGN continuum variability, and it is argued 
here that such variable emission also readily explains a 
wide variety of hitherto puzzling observations of AGNs. 
Things explained include: 

1. the different lags measured for different events in 
reverberation-mapping campaigns, 

2. the unusually asymmetric line profiles seen in some 
AGNs, 

3. line-profile variability (including sharp structures 
seen in difference and rms spectra), 

4. the strange line flux and continuum correlations or 
non-correlations often seen over remarkably small 
velocity intervals, 

5. the sometimes rapidly changing and conflicting 
kinematic signatures seen in velocity-resolved re- 
verberation mapping, 

6. why variable bumps in line profiles are more pro- 
nounced in H/3 than Ha, 



7. velocity-dependent line polarization 

8. line-polarization variability 

An important related conclusion is that a single BLR 
geomet ry with the structure advocated by iGaskell et al.l 
(2007b]), and with kinematics consistent with this, suc- 
cessfully explains a wide variety of BLR observations 
once the presence of off-axis variability is recognized. 
AGNs with very different line profiles are not fundamen- 
tally different. In particular, AGNs showing disk-like line 
profiles are not fundamentally different from A GNs show- 
ing classical logarithmic profiles. 

It is clear that orientation has a significant effect on line 
profiles. This will introduce systematic errors in black 
hole mass estimates. 

While more exotic explanations of the various phenom- 
ena considered here are not necessarily ruled out (e.g., 
binary black holes, recoiling black holes, or warped ac- 
cretion disks) , the off-axis variability model renders them 
unnecessary. Off-axis variability also makes it very hard 
to detect such phenomena from observations of broad 
lines. 

The complications created by off-axis variability need 
to be taken into account when interpreting reverberation- 
mapping results and planning monitoring campaigns. 
These complications severely limit what can be learned 
from reverberation mapping. In particular the goal of 
detecting fine structure in BLRs seems unattainable. 
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